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Abstract: 

To match the precision of present and future measurements of jy-boson produc- 
tion at hadron colliders electroweak radiative corrections must be included in the 
theory predictions. In this paper we consider their effect on the transverse momen- 
tum (pr) distribution of W bosons, with emphasis on large pp- We evaluate the full 
00 ■ electroweak O(a) corrections to the processes pp — *> W + jet and pp — > W + jet in- 

cluding virtual and real photonic contributions. We present the explicit expressions 
in analytical form for the virtual corrections and provide results for the real correc- 
tions, discussing in detail the treatment of soft and collinear singularities. We also 
provide compact approximate expressions which are valid in the high-energy region, 
where the electroweak corrections are strongly enhanced by logarithms of s/M^r. 
These expressions describe the complete asymptotic behaviour at one loop as well 
as the leading and next-to-leading logarithms at two loops. Numerical results are 
presented for proton-proton collisions at 14 TeV and proton-antiproton collisions at 
2 TeV. The corrections are negative and their size increases with pp- At the LHC, 
where transverse momenta of 2 TeV or more can be reached, the one- and two-loop 
corrections amount up to —40% and +10%, respectively, and will be important for 
a precise analysis of W production. At the Tevatron, transverse momenta up to 
300 GeV are within reach. In this case the electroweak corrections amount up to 
— 10% and are thus larger than the expected statistical error. 



1 Introduction 



After the startup of the Large Hadron Collider (LHC) hard scattering reactions will 
be explored with high event rates and momentum transfers up to several TeV. In 
order to identify new phenomena in this region, the predictions of the Standard 
Model have to be understood with adequate precision. 

The study of gauge-boson production has been among the primary goals of 
hadron colliders, starting with the discovery of the W and Z bosons more than two 
decades ago pQ. The investigation of the production dynamics, strictly predicted 
by the electroweak theory, constitutes one of the important tests of the Standard 
Model. Differential distributions of gauge bosons, in rapidity as well as in transverse 
momentum (pr), have always been the subject of theoretical and experimental stud- 
ies. This allows to search for and set limits on anomalous gauge-boson couplings, 
measure the parton distribution functions and, if understood sufficiently well, use 
these reactions to calibrate the luminosity. For gauge-boson production at large pr 
the final state of the leading-order process consists of an electroweak gauge boson 
plus one recoiling jet. Being, in leading order, proportional to the strong coupling 
constant, these reactions could also lead to a determination of as in the TeV region. 

The high center-of-mass energy at the LHC in combination with its enormous 
luminosity will allow to produce gauge bosons with transverse momenta up to 2 TeV 
or even beyond. In this kinematic region the electroweak corrections are strongly en- 
hanced, with the dominant terms in L-loop approximation being leading logarithms 
(LL) of the form a L log 2L (s/M^), next-to-leading logarithms (NLL) of the form 
a L log 2L_1 (s/M^/), and so on. These corrections, also known as electroweak Sudakov 
logarithms, may well amount to several tens of percent [21E1SE1EIEIIEJE1 EH E] ■ 
(A recent survey of the literature on electroweak Sudakov logarithms can be found 
in Ref. [12J.) Specifically, the electroweak corrections to the pr-distribution of 
photons and Z bosons at hadron colliders were studied in Refs. [7J El [91 \10\ . In 
Refs. [HI El CD], it was found that at transverse momenta of (9(1 TeV) the dominant 
two-loop contributions to these reactions amount to several percent and must be 
included to match the precision of the LHC experiments. This is quite different 
from the production of on-shell gauge bosons with small transverse momenta [13J, 
where the electroweak corrections are not enhanced by Sudakov logarithms. With 
this motivation in mind we study the electroweak corrections to hadronic production 
of W bosons in association with a jet, pp — > Wj, at large pp- 

As a consequence of the non-vanishing W charge, QED corrections cannot be 
separated from the purely weak ones and will thus be included in our analysis. Thus, 
in comparison with Z-boson production, several new aspects arise. Real photon 
emission must be included to cancel the infrared divergencies from virtual photonic 
corrections. Collinear singularities, a consequence of radiation from massless quarks, 
must be isolated and absorbed in the parton distribution functions (PDFs) in the 
case of initial-state radiation. We regularize soft and collinear singularities in two 
different schemes: using small quark and photon masses which are set to zero at the 
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end of the calculation and, alternatively, dimensional regularization. In events with 
real radiation, the pt of the W boson is balanced both by the p<r of the recoiling 
parton (quark or gluon) and the photon. Configurations involving a small-pr parton 
and a hard photon are better described as Wj final states. We thus define the Wj 
cross section imposing a lower limit on the jet transverse momentum, which is chosen 
independent of the W-boson pr- in order to avoid final-state collinear singularities, 
we recombine collinear photon-quark final states. 

The virtual EW corrections to Wj production are formally connected with the 
real emission of W and Z bosons, which leads to WVj final states with V = W, Z. 
Both contributions are of 0(a 2 as)- If integrated over the full phase space, the real 
emission of gauge bosons produces large Sudakov logarithms that partially cancel 
those resulting from virtual gauge bosons. However, in exclusive measurements of 
pp — > Wj, the available phase space for gauge boson emission is strongly suppressed 
by the experimental cuts. We thus expect that real emission provides relatively 
small contributions while the bulk of electroweak effects originates from virtual 
corrections. In fact, for pp — > Zj it was shown that, in presence of realistic (and 
relatively less exclusive) experimental cuts, the contribution of real emission is about 
five times smaller than the virtual corrections [11] . Moreover, real emission can be 
further reduced with a veto on additional jets, which suppresses multiple-jet events 
resulting from the hadronic decay of the radiated gauge bosons. Therefore we will 
restrict ourselves to the investigation of virtual electroweak corrections (and photon 
bremsstrahlung) . The real emission of W and Z bosons can be non- negligible and 
certainly deserves further detailed studies, however we do not expect a dramatic 
impact on our results. 

The partonic reactions qq' — > W ± g( r y), q'g — > W ± q (7) and qg — > W ± q' (7) 
with q = u,d, s, c, b are considered. All of them are, however, trivially related by 
CP- and crossing-symmetry relations. Quark-mass effects are neglected throughout, 
which allows to incorporate the effect of quark mixing through a simple redefini- 
tion of parton distribution functions (see Sect. 12.11) . Our conventions for couplings, 
kinematics and two- as well as three-body phase space are also collected in Sect. 2. 
The calculation of the virtual corrections is described in Sect. El We present an- 
alytic expressions for the one-loop amplitude, specify the counterterms in the 
renormalization scheme and isolate the infrared singularities. The high-energy limit 
is studied in detail in Sect.Hl The analytic one- loop result is investigated in the 
limit s ^> Myy, keeping quadratic and linear logarithms as well as constant terms. 
These results are compared to those derived in the NLL approximation [T3] . In view 
of their numerical importance we also derive the dominant (NLL) two-loop terms, 
using the formalism of Refs. [T5l [T6] . The calculation of the real corrections is per- 
formed using the dipole subtraction formalism [T7J, UHl EE], as discussed in Sect. [5] 
The checks which we carried out in order to ensure the correctness of the results are 
described in Sect. El 

The numerical results are presented in Sect. [71 After convolution with parton dis- 
tribution functions, we obtain radiatively corrected predictions for pp-distributions 
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of W bosons at the LHC and the Tevatron. The quality of the one-loop NLL and 
NNLL approximations is investigated and the size of the dominant two-loop terms 
is compared with the expected statistical precision of the experiments. Concerning 
perturbative QCD, our predictions are based on the lowest order. To obtain real- 
istic absolute cross sections, higher-order QCD corrections |20j must be included. 
However, the relative rates for W + , W~ and Z production are expected to be more 
stable against QCD effects. Therefore, the impact of the electroweak corrections on 
these ratios is presented in Sect. [7J Our conclusions and a brief summary can be 
found in Sect. [HJ Explicit analytic results are collected in the Appendices. 

A short description of the method of calculation and the main results for LHC 
have been given in Ref. [2T] . After completion of this work, Hollik, Kasprzik and 
Kniehl [22J reported results on hadronic H^-boson production at large Pt qualita- 
tively similar to those of Ref. [21] and the present paper. 

2 Definitions and conventions 
2.1 Hadronic cross section 

The pr-distribution of W bosons in the reaction hiti2 —>■ W ± j( , y) is given by 



— ] = Y] / dxi / dx 2 9{xix 2 -f iairi )f hlA (xi 1 ii 2 )f h2 ^(x2^ 2 ) . 



where f m i n depends on the kinematic configuration of the final state and is specified 
at the end of Sect. 12.21 The indices a, b denote initial-state partons and fh lt a(x, /i 2 ), 
fh 2 ,b{ x ) A* 2 ) are the corresponding parton distribution functions (PDFs). (j ab ^ wak (i) 
is the partonic cross section for the subprocess ab — > W a k( , y). The sum in ([T]) runs 
over all a, b, k combinations corresponding to the subprocesses 

d n u m -> W + g('j), u m d n -> W + g('j), gu m -> W+dJfl), 

u m g -> W + d n {-i), d n g iy + w m (7), gd n -> W + u m {"i), (2) 

for W + production, and similarly for W~ production. 

The dependence of the partonic cross sections on the family indices m, n amounts 
to an overall factor |V^ m d n | 2 . This factor can be easily absorbed by redefining the 
parton distribution functions as 

3 3 

fh,d m = I^Umdnl fh,dn1 fh,dm = l^mrfnl fh,d n 1 

71=1 n=l 
fh,u m = fh,Umi fh,u m = fh,u m i fh,g = fh,g- (3) 

The hadronic cross section (pQ) can be computed using the trivial CKM matrix 
Vuidj = $ij and the redefined PDFs (j3J). Since we do not consider initial or final 
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states involving (anti-)top quarks, only 
families (m = 1, 2) have to be included, 
with q = u,d, c, s) automatically include 
initial- and final-state bottom quarks. 



the contributions of the first two quark 
The corresponding redefined PDFs {fh,q 
the (small) contributions associated with 



2.2 Kinematics 

For the 2^2 subprocess ab — > W c k the Mandelstam variables are defined in the 
standard way, 

s = (p a + Pb) 2 , t=(p a -p w ) 2 , u={p b -p w f. (4) 

The momenta p a , p b , Pk of the partons are assumed to be massless, whereas ply = 
M^y. In terms of Xi,x 2 ,pr and the collider energy a/s we have 

s = x\X2S, i= — — (1 — cos6), u= — — (l + cos6), (5) 



with cos 6 = J 1 — <ip\s/(s — My/) 2 corresponding to the cosine of the angle between 
the momenta p a and pw in the partonic center-of-mass frame. 

The pr-distribution for the unpolarized partonic subprocess ab — > W a k reads 



dd ab ~ 



■Wk 



dp T 



ab 



J d$ 2 £|Ar^T^o, 2 ($ 2 ), (6) 



where = \ J2 po \ J2 co \ involves the sum over polarization and color as well as the 
average factor 1/4 for initial-state polarization. The factor M a b is given by 

sr _ (^) 4 m 

where N qq > = N qq , = N*, N gq = N qg = N qg = N gq = N C (N 2 - 1), with iV c = 3, 
account for the initial-state colour average. The phase-space measure d<3>2 is given 
by 

d 3 p w d 3 p k 4 
d$2 = (2vr)32^ ( 27 r)32rf 5 "^ "P*) , ^ 

while the function Fq 2 defines the observable of interest, i.e. the W-boson p-^- 
distribution in presence of a cut on the transverse momentum of the jet, 

FoA^) = S(p T -p T>w )9(p TJ -p$) . (9) 

In the 2-particle phase space the jet is identified with the parton k and momentum 
conservation implies ptj = pT,k = Pt, w- In practice, since we always consider the 
Px-distribution in the region p T > > Ptji the cut on p T j in (J5j) is irrelevant. 
The phase-space integral in yields two contributions originating from kinematic 
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configurations in the forward and backward hemispheres with opposite values of 
cos 9 in the center-of-mass frame, 



J-ab^Wk J-ab-^Wk J-ab^Wk 
u " . <l(T Iw(l | QO bkwd 



dpT dpx dp 

with 



A-ab-^Wk „ j-ab^Wk j-ab^Wk 

(l "i\vd /'I ^VJ"'' (l "l)k\v<l (l "l\vd 



dp T 8nN ab s\t-u\ ' dp T dpq 



.(11) 



For the 2^3 subprocess a& — > W a k'y we define the following five independent 
invariants 

s = (p a + Pb) 2 , i = (p a - p w ) 2 , u= (p b - p w ) 2 , 

i' = (Pa-p 7 ) 2 , u'=(p b - Pl )\ (12) 

and the four dependent invariants 

s' = (p k + P 7 ) 2 = s + i + u- M^, s" = (p w + p k f = s + t' + u', 
i" = (Pa - Pk) 2 = M 2 v -s-t- t', u" = (p b - p k f = M$y-s-u-u' 



(13) 



The pr-distribution for this subprocess reads 

^ab~*W"k^ 

dp T 

where 



Nab J d$ 3 Y}M ab ^ w ^\ 2 Fo, 3 ($3), (14) 



,, d 3 p w d 3 p k d 3 p 7 4 

d$3 = (2^2^ ( 27 r)32p0 5 (P-+P6 -p* . (15) 

In the 3-particle phase space, the W-boson px-distribution in Wj production is 
defined by the observable function 

FoA^s) = 5(PT-PT,w)9( P T,j-ptj) ■ (16) 

The cut on the jet transverse momentum rejects events where the W^-boson px is 
balanced by an isolated photon plus a parton with small transverse momentum. This 
observable is thus free from singularities associated with soft and collinear quarks or 
gluons. When applying the cut on the jet momentum in the 3-particle phase space, 
care must be taken that the definition of the jet px is collinear-safe. In general the jet 
cannot be identified with the parton k, since in presence of collinear photon radiation 
the transverse momentum of a charged parton is not a collinear-safe quantity. Thus 
we identify the jet with the parton k only if k is a quark well separated from the 
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photon or a gluon. Otherwise, i.e. for collinear quark-photon configurations, the 
recombined momentum of the quark and photon is taken as momentum of the jet. 
In practice, we define the separation variable 



R^l) = ^(r lq - V ,f + (<p q -<P 1 )\ (17) 

where rji is the pseudo-rapidity and fa is the azimuthal angle of a particle i. If 
R(<1j1) < -Rsep, then the photon and quark momenta are recombined by simple 



four- vector addition into an effective momentum pj and then px, j = y (PT,q + Pt, 7 ) 2 , 
otherwise p^j = pr,q- We note that, in the collinear region, lowest-order kinematics 
implies ptj = px.q + Pt, 7 = Pt,w > Ptj- This means that the recombination 
procedure effectively removes the cut on p T ^ inside the collinear cone R(q, 7) < -R sep - 
For instance the recombined gq' — > W a q r y cross section is given by 

°&r W °' n = I da 9 *'-*™^ + I 6(p T q -p^)da 9q '^ w ^\ (18) 

J R(q,-y)<R scp JR(q,-y)>R scp 

In contrast, for the case of final-state gluons, we do not perform photon-gluon re- 
combination and the cut on p T> g is imposed in the entire phase space. 

This procedure has the advantage to avoid both collinear-photon and soft-gluon 
singularities. However it implies a different treatment of quark and gluon final states 
and can thus be regarded as an arbitrary cut-off prescription for the final-state 
collinear singularity. Moreover, the recombined cross section ({TBI has a logarithmic 
dependence on the cut-off parameter -R sep . These aspects are discussed in detail 
in Appendix [XI There we compare the recombination procedure with a realistic 
experimental definition of exclusive pp — > Wj production, where final-state quarks 
are subject to the same cut as final state gluons (pT, q > Ptj) within the entire phase 
space. Describing the exclusive gq' — > W a q^ cross section, 

a 9 iZ W ^ = J 9(p T , q - p™) da 9 ^ w ^, (19) 

by means of quark fragmentation functions, we find that the quantitative difference 
between the two definitions (lisp and (fl^j) amounts to less than two permille. More- 
over, we show that the recombined cross section is extremely stable with respect 
to variations of the parameter i? S ep- This means that the recombination procedure 
used in our calculation provides a very good description of exclusive pp — > Wj 
production. 

Another treatment of the singularities, which does not require recombination 
and treats quark- and gluon-induced jets uniformly, has been proposed in Ref. [22J. 
There, contributions from Wj production and W^ production to a more inclusive 
observable, i.e. high-px W production, are both calculated. All soft and collinear 
singularities in the final state cancel in the approach of Ref. [22] as a result of the 
more inclusive observable definition than associated production of the W boson 
together with a jet, considered in this work. The comparison of our results with 
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those of Ref. [22] seems to indicate that these differences in the jet definitions have 
a quite small impact on the size of the electroweak corrections. 

The quantity f min in (pQ) is related to the minimum partonic energy that is needed 
to produce final states with p T j > p^j and Pt,w > Pt'wo 

sf miQ = (p™ + ^(p^ w ) 2 + M&) 2 ■ (20) 

When we evaluate the 2^2 contributions to the hadronic cross section (jrj, after 
analytic integration of the phase space in (jSJ), we can set p™™ = p™w — Pt in (1211 . 



2.3 Crossing symmetries 

The unpolarized squared matrix elements for the 2 — * 2 processes in are related 
by the crossing-symmetry relations 



^\ i M 9q '^ w " q \ 2 = - ^2\M g ~ q '^ w ' T9 \ 2 
^2\M q9 ^ w<Tq '\ 2 = - ^\ l M qq '^ w ' J9 \ 2 



s<->t 



(21) 



t^rU 



Moreover, due to CP symmetry, the unpolarized partonic cross section for the pro- 
duction of positively and negatively charged W bosons are related by 



[22) 



Eqs. ( 1211) and ( 1221) permit to relate the six processes for W + production in (jzj) and 
the six charge conjugate ones to a single process. Hence the explicit computation 
of the unpolarized squared matrix element needs to be performed only once. In the 
following we will present explicit results for the process qq' — > W a g. 

Similarly, for the unpolarized squared matrix elements for the 2^3 processes 
in (121) we have 



gq'^Wqj\2 



q-q'->W°9J\2 



Y^\M g ~ 9 ^ w ' Tq ~'' y \ 2 = - Y^\M 9 ~ q '~* w ' T9 ' y \ 2 



{$^>u" ,t++s" ,t' *->$'} 
{s^i",u^s",u' <->s'} 



and 



It is thus enough to perform calculations only for the qq' — > W cr g , ~f subprocess. 



(23) 
(24) 
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2.4 Couplings and Born matrix element 

For gauge couplings we adopt the conventions of Ref. [23]. With this notation the 
gqq vertex and the Vq'q vertices with V = A, Z, W ± read 

q / q' / -\-'<-L 

where uj\ are the chiral projectors 



rV f 
"<A<A' 



(25) 



1, s 1, 

-(l+ 75 ), WL = -fl 



2 , - 2 75), (26) 

t a are the Gell-Mann matrices and I v are matrices in the weak isospin space. For 
diagonal matrices such as I A and I z we write I v , = 8 aa 'lY . In terms of the weak 



isospin and the weak hypercharge Y qx we have 



<a 



Q 



Y 

_rp3 _ 1 Qx 

«* 2 ' 



,. ■) '/.\ 'A •) {"^) 

with the shorthands c w = cos^w and s w = sin#w for the weak mixing angle #w 
The eigenvalues of isospin, hypercharge and SU(2) Casimir operators for left-handed 
fermions are 



T 



-Ti 



I Y 

2' 



3, c F 



4' Ca 



(2f 



The only non-vanishing components of the generators associated with W bosons are 

1 



jW+ 
u^d L 



jW~ 



(29) 



The triple gauge-bosons vertices read 




£ V a V b V c ^gVW** ( 

+ 9 w (fc3-^n (30) 

where the totally anti-symmetric tensor e lV2Vs is defined through the commutation 
relations 



j Vl rV 2 



and has components e zw+w = — ic w and e / "' " = is w . 

To lowest order in a and as, the unpolarized squared matrix element for the 
qq' — > W a g process reads 

Y i \m$-* w " b \ 2 



£ VlV2Vsjv 3 

V 3 =A,Z,W± 



(31) 



-AW+W- 



H 2 aa s (N 2 -l) (i^r) 



2 t 2 + u 2 + 2M^s 



tii 



(32) 



where a = e 2 / (Att) and a$ = 9$/ (47r) are the electromagnetic and the strong coupling 
constants. 
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MAAAA/ 



WWW 



(tl) 



(t2) 



Figure 1: Tree-level Feynman diagrams for the process qq' — > W a g. 

3 Virtual corrections 

In this section we present the virtual electroweak corrections to the qq' — > W a g 
process. The algebraic reduction to gauge-coupling structures, standard matrix 
elements and one-loop scalar integrals is described in Sect. 13.21 The renormalization 
of ultraviolet divergences and the subtraction of infrared singularities originating 
from soft and collinear virtual photons are discussed in Sect. 13.31 and Sect. 13.41 
respectively. In Sect. 13.51 we summarize the one-loop result for the unpolarized 
squared matrix element. 



3.1 Preliminaries 

As discussed in the previous section, the twelve different processes relevant for Wj 
production are related by CP and crossing symmetries. It is thus sufficient to con- 
sider only one of these processes. In the following we derive the one-loop corrections 
for the qq' — > W cr g process. The matrix element 

Mf' ~* W ' B = M^ w ' a + 5Ml q '^ 9 (33) 

is expressed as a function of the Mandelstam invariants 

S=(Pq+Pq') 2 , t=(Pq-Pw) 2 , « = {Pq> ~ Vwf ■ (34) 

The Born contribution Ai^ 9 ~* w 9 results from the t- and u-channel diagrams of 
Fig. [H The loop and counterterm diagrams contributing to the corrections, 

5M f'^ 9 = 5M f^; 9 + 6Mt& w '; (35) 

are depicted in Fig. [2] and Fig. [31 respectively. 

The quarks that are present in the loop diagrams of Fig. [2] are treated as massless, 
and the regularization of the collinear singularities that arise in this limit is discussed 
in Sect. 13.41 The only quark-mass effects that we take into account are the m t - terms 
that enter the counterterms through gauge-boson self-energies. 

Our calculation has been performed at the matrix-element level and provides full 
control over polarization effects. However, at this level, the analytical expressions 
are too large to be published. Explicit results will thus be presented only for the 
unpolarized squared matrix element. 
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(bl) 



(b2) 



(b3) 



Figure 2: One-loop Feynman diagrams for the process qq' — > W a g. The diagrams 
si, s2, vl and v2 receive contributions from neutral and charged gauge bosons, 

V = A, Z, W . The diagrams v3, v4, bl and b2 involve only neutral gauge bosons, 

V = A, Z. The remaining diagrams, v5, v6 and b3 involve two contributions with 
one charged and one neutral gauge boson: (Vi, V2) = (V, and (W ± , V) with 

V = A,Z. 




SLSLSJLQSls 




(cl) 



(c2) 



iWWVV 



(c3) 



www^ 



(c4) 



®AAAA/V 



-kSlSLSSLSJls 

c5) 



■#/vwv 
(c6) 



Figure 3: Counterterm diagrams for the process qq' — > W a g. 
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3.2 Algebraic reduction 

The matrix element ( |33i) has the general form 

M 9f^w a = ieg^vip^M^uMP^^iPwKiPs). (36) 

Since we neglect quark masses, jW^'^ consists of terms involving an odd number of 
matrices 7 P with p — 0, . . . , 3. The 7 5 -terms are isolated in the chiral projector 
defined in (I26p . The polarization dependence of the quark spinors and gauge-boson 
polarization vectors is implicitly understood. In analogy to (l33l) and (l35l) we write 

MY V = MY" + SM^, bMY v = SM^Zps + (37) 

Following the approach adopted in Ref. |9j, we isolate the SU(2)xU(l) couplings 
that appear in the Feynman diagrams and reduce the one-loop amplitude to a sum 
of contributions associated with independent coupling structures. As we will see, 
besides an abelian and a non-abelian contribution that are related to the ones found 
for Zj production [9], for Wj production we have two additional coupling structures. 

The coupling structure of the Born amplitude is trivial and consists simply of 
the QhQij component of the SU(2) generator, 

M =I qhqi S =- 7 =-, S = ? + . (38) 

The contribution of the loop diagrams of Fig. [2] can be written as 



47T [ 



y=A,z,w± 



(l w -°I v f) Df(A£) + (l v I v I w ~ a ) Di?{M v ) 



+ E 



y=A,z 



(l v I w -'l v \ D%"(My) + ±e wavw -' (l v I l 

V / m.ai .<!„. V 



x D%"(My, M 2 W ) + J- £ vw°w-° fjw-jv} D^{M 2 W1 My 



(39) 



In the following, treating the electroweak gauge couplings as isospin matrices and 
using group-theoretical identities (see App. B of Ref. [23]), we express the above 
amplitude in terms of the eigenvalues of isospin, hypercharge and SU(2) Casimir 
operators for left-handed fermions (1281) . 

The tensors D^ 1 '{My) and D 1 ^ {My) in (1391) describe the contributions of the 
diagrams si, vl and s2, v2, respectively. These diagrams may involve charged or 
neutral virtual bosons. In the former case (V = W^, the corresponding couplings 
reacQ 

E i w ~ a i WP i WP = e i wp i w ~ p i w ~ a = ^ F ~ (-^ ,3 ) 2 jw-° ^q-j 

p=± p=± s w 



1 The following identities have to be understood as matrix identities, where the <7l9l indices of 
the SU(2) generators are implicitly understood. 
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In the latter case (V = A, Z) the coupling factors read 

(T 3 )2 



rW- 



, I V I V 



SU(2) 
VV 



jVjVjW- 



SU(2) 
VV 



pay 



Here 



? SU(2) 
'VV 



{Uvw 3 ) 



Xx 



+ X V T Y + 5 
- X V T 3 Y 4 S 

Uvw 3 UvB 

s w c w 



U(l) 


Y 2 1 


VV 


4c 2 


U(l) 




VV 


4c 2 



rW~ 



fU(l) 

Vv 



ysj 



(41) 



(42) 



where U is the electroweak mixing matrix. For V 

X A = -1, 5^ = c^, and <5§^ (2) = c^, X z = 1, 5^ = s 2 „. We note that 



A, Z we have S^ 2 ^ 



8 

V=A,Z 



SU(2) 
VV 



E * 

V=A,Z 



U(l) 

yy 



E ^ 



(43) 



V=A,Z 



The tensor D^ U (M V ) in ( |39l) corresponds to the diagrams v3, v4, bl and b2. These 
diagrams receive contributions from neutral virtual gauge bosons only (V = A, Z). 
For the corresponding couplings we have 



V tW~ 



ri 



rV 



SU(2) 
VV 



Cf — Ca/2 



+ 5 



u(i) Y 2 



(44) 



Finally, D% v (M Vl ,M V2 ) represents the dia grams v5, v6 and b3. These diagrams 
involve a neutral gauge boson (V = A, Z) and a W boson. The coupling factors 
yield 



avvw- 



V tW~ 



I V I 



4 u v ( 2 )^4-XvT 3 F 



1 ^WW-" jW-" jV 



A SU(2) W A . y T 3 Y 



4s 2 
C A 



jiy- 



(45) 



Using the above identities we express the one-loop amplitude ( I3U1) for lU-boson 
production in a form that is analogous to the one adopted in Refs. 0, [10] to describe 
the production of neutral gauge bosons. To this end we defin^] 



Dl U (My) + D^(My) 



D^(My), 



V) 



DT{My, M 2 W ) + ^(4, My) - D^(My) 



5A^(M V ) = D^(M 2 V ) + D?(M\ 
SAtUMy) = DT(M V} Ml) - DT(M 2 W} My) + D? (M 2 ) - D^\M 2 ). (46) 



2 In our notation we emphasize the dependence of the form factors 5A^ \ on My, whereas the 
dependence on the external momenta as well as the Mjy-dependence (for 5A^ V N and SA^y) is 
implicitly understood. 
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The tensor 8Ai\(M v ) is identical to the abelian tensor defined in Ref. [9], and 
5Ai u N (My) is equal to the non-abelian tensor of Ref. [9] for M v = Myy- The re- 
maining two tensors, 5Ai\(M v ) and 6Aiy(M v ), are new. Using (T4T)]) - (14"6]) we can 
write the one-loop amplitude 



as 



5M\ 



loops 



a 



4tt\/2.< 



E 



w KV=A,Z 



6 



SU(2) 

vv 



c 



_ 5 SU(2) Cf 



'VV 



w / 

2 SA^(M V ) 



p^5Ar x (My)+s v T )CA 



- X v T*Y q JA^(M } 



+ 



- (7* 



9L> 



2s 



<Uf x (M, 



w j 



(47) 



This amplitude has been reduced algebraically using the Dirac equation, the identity 
p^Efj.ip) = for gauge-boson polarization vectors and Dirac algebra. Moreover, 
tensor loop integrals have been reduced to scalar ones by means of the Passarino- 
Veltman technique [21] . The result has been expressed in the form 



10 



5az(m v ) = Y.Y.n\My)sruM, 



v), 



(4* 



1=1 j 



for I=A,N,X,Y. The quantities ^(My) are rational functions of Mandelstam in- 
variants and masses. Explicit expressions for the tensors Sf v and the scalar loop 
integrals Jj(M v ) are provided in Appendix IB1 and Appendix O 



3.3 Renormalization 

While the tensors 5A X x and 5Ai Y are ultraviolet finite, the abelian and the non- 
abelian tensors give rise to the ultraviolet singularities 



5ATM 



uv 



Attv S, 



UV^o ) 



uv 



2A UV S 



(49) 



where Sq U is the tensor structure of the Born amplitude ( 1381) . and 



'47r/i 2 V r(l+e) 



Ml 



- - 7e + ln(4vr) + In ' '" 



Ml 



+ 0(e) (50) 



in D = 4 — 2e dimensions. These singularities are cancelled by the counterterm 
diagrams depicted in Fig. [3] and the results are independent of the scale \x of dimen- 
sional regularization. The counterterms that are responsible for the contributions 
of diagrams cl, c2, c3 and c4 read 



3d 



^^0o~5757P 



-i^VW^L- (51) 
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Since there is no 0{a) contribution to the renormalization of the strong coupling con- 
stant gs, these counterterms depend only on the wave-function renormalization con- 
stants for left-handed quarks, 5Z qL . Their combined contribution to the qq' — »■ W a g 
process, i.e. the sum of the diagrams cl, c2, c3 and c4, vanishes. The renormaliza- 
tion of the qq 1 — > W a g process is thus provided by the diagrams c5 and c6, which 
originate from the Wqq' counterterm, 



with 




w~ a 



5C A + 6C* 



(52) 



5C 1 



(sz„ 



"L 



sz, 



5C r 




and yields 



5M\^ T = (5C A + 5C N ) M 







(53) 



(54) 



The wave-function renormalization constants of massless left-handed quarks and 
on-shell W bosons are related to the corresponding self-energies by 



5Z qL = -Re (E 9 ' L (0) 



5Z 



w 



-Re 



dp 2 



(55) 



and have been evaluated using the explicit results of Ref. 

For the definition and the renormalization of the SU(2) coupling constant, 



9l 



4ira 5g\ 5a 



5s 2 



9 ' 2 9 ' 

< 92 « si 



(56) 



we adopt the G^-scheme, where the electromagnetic coupling constant a is expressed 
in terms of the Fermi constant G M , and the weak mixing angle is related to the on- 
shell masses Mz, My/ of the gauge bosons, 



a 



l-ci 



IT 



1 - M^/Mf . 



(57) 



The counterterm 5a /a in the G M -scheme can be derived from the on-shell counter- 
term 5a(0)/a(0) for the fine-structure constant in the Thompson limit. Using the 
one-loop relation a = a(0) [1 + Ar] and requiring a + 5a = a(0) + 5a(0) we have 

5a 5a(0) 
a 



a(0) 



Ar. 



(5* 



Combining the relations 
l26l [271. we obtain 



and using the explicit one-loop expression for Ar 



5gl 



-Af = Re 

92 



£^(M^)-£^(0) 



M 2 W 



a 



A 



uv 




7-124 X JM 2 W 



2 s 2 



(59) 
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The above conterterms yield the ultraviolet singularities 



5C A 
SC N 



uv 



uv 



a 



2-ksIj 



A 



uv 



(60) 



Using ( )43l) one can easily verify that these singularities cancel those resulting from 
the loop diagrams [see ( 14"7T) and ( l4~9i) ]. 



3.4 Soft and collinear singularities 

Loop diagrams and wave-function renormalization constants involve singularities 
originating from soft and collinear virtual photons (for brevity denoted in the fol- 
lowing as IR singularities) . In order to isolate these singularities and check that they 
are cancelled by corresponding ones originating from real photon bremsstrahlung, 
we split the wave-function renormalization constants and the photon contributions 
to (H71) . i.e. the terms (L4^j(M|), in IR-singular (IR) and IR-finite (fin) parts: 

5Z qL = + 62%, 
6Z W = 6Z% + 5Z$, 
5AZ(M 2 A ) = 6Af/ lu + 8A\y v . (61) 

The singular parts depend on the scheme adopted to regularize IR singularities. 
The remaining parts are scheme-independent and free from IR singularities, but 
can contain ultraviolet poles. For the regularization of IR singularities we use, 
alternatively, two different schemes: 

• In the first scheme, which we denote as mass-regularization scheme (MR), we 
use infinitesimal quark masses m and a photon-mass regulator, Ma = A with 
< A <C m. Since the quark- mass dependence disappears in the final result, 
we perform the computation using the same mass m for all quarks. To denote 
quantities evaluated in this scheme we use the label MR; 

• In the second scheme we perform the calculation using massless fermions and 
photons, Ma = m = 0, and we evaluate IR singularities in dimensional reg- 
ularization (DR). To denote quantities evaluated in this scheme we use the 
label DR. 



The singular parts of the wave-function renormalization constants read 



7U , ^ a ( Anfx 2 , m 



^w = T-\^r) r(i + e)C (62) 



47T \M^ 
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with 



ij, MR 



2 In 



A 2 



in the MR scheme and 



LlR 

<7,DR 



1 



-4, 



/IR 

W.DR 



W.MR 



-2 In 



A 2 



(63) 



(64) 



e ' e 

in the DR scheme. The splitting of the loop contributions (L4^j(M|) into IR-singular 
and IR- finite parts is performed at the level of the scalar loop integrals Ji(M\): 

mi) = ji r + jf n . (65) 



Explicit expression for the IR-singular and IR-finite parts of individual loop integrals 
are presented in Appendix [D] Combining all singular contributions Jj R we obtain 

m,»u _ ( 4?r/x 2X 



M w 



IR Q/J-f 
5 



(66) 



i.e. the IR singularities factorizqj with respect to the Born amplitude (]38|) . The 
IR-singular part of the renormalized amplitude can be expressed in terms of the 
electromagnetic charges of the external particles as 



M 2 W 



r (1 + e) 



Q q Q q >fl R + vQj2 K -°Q q >f: 



IK 



cIR 



MY\ (67) 



/■IR 
' Y 



where a = ±1 is the charge of the W boson and 

A IR = -/I R - hf, 

f lR flR flR , f rlR _ f 

J2 — J A J N 2 \ X 

rlR _ rlR f IR , ^ f flR , flR >TR l 
i3 — ~ J A ~ J N + 2 V X ^ Y ~ 9 _ 

In the MR scheme we obtain 
/ A 2 



IR 



h lti — h l 



(6£ 



fi 



IR 

,MR 



— 2 In 



— s 



M w. 



In — - W 



flR 
J'. 



2,MR 



flR 
i3,MR 



hi 



A 2 



in 



In 

2 \ 



mf 



m 



+ 2 In 



m 



+ 3 In 



2 In ^1 

A 2 \ 



-In 2 



+ 21n| — 



m 



M w , 



flR 

J 2, MR 



t— m 



(69) 



3 To be precise, the tensors Mj" x {M\) and <5^4^(M^) contain also non-factorizable IR diver- 
gences. However these non-factorizable singularities are related by 



2M™ 1 ~ fact <H^), 



and due to the identity Cf — (T ( 



3 \2 _ 



the X- and N-terms in (jTf|) . they cancel. 



Ca/4, which relates the coupling structures associated with 
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and in the DR scheme 



Jl.DR 



J2,DR 
i3,DR 



1 

£ 
1 

e 



2 In 



My/. 



2 In 1 



My/, 



-4, 

5" 
2 



IR 

2,DR 



t— >M 



(70) 



The splitting ( 1611) has been performed in such a way that in the high-energy limit 
(s, \t\, \u\ My/) the IR-finite part of the amplitude has the same logarithmic be- 
haviour as the virtual corrections regularized by a photon mass = My/. Indeed 
the IR-singular parts fj R correspond exactly to the contribution called purely elec- 
tromagnetic in Ref. [H]. This implies that, up to terms that are not logarithmically 
enhanced at high energies, the IR-finite part of the corrections corresponds to the 
symmetric electroweak contribution of Ref. [H] , which is constructed by setting the 
photon mass equal to My/. This property is evident in the asymptotic high-energy 
expressions ( 1801) for the IR-finite part of the diagrams involving virtual photons. 



3.5 Result 



Let us summarize our result for the unpolarized squared matrix element for the 
qq' — > W a g process. To 0(a 2 as), 



+ 2Re 



qq'^Wg 



. (71) 



Using (I3T)]) and summing over the polarizations we can express the interference term 

as 



2Re 



qq'^Wg 



< Re [Tr (^M^^M^' 



2w 2 aa s (N 2 - 1) 



9w' g^i' 



Pw 



(72) 



where Ai = 7°A^^7°. Combining the contributions of the bare one- loop diagrams 
( 1471) and the counterterms (154"1) yields 

Y} M f' -+ W ° 9 \ 2 = [l + 2 Re (5C A + Y}Mf' ^ W ° 9 \ 2 + ^^(N 2 - 1) 



xRe £ 

lv=A,Z 



' rSU(2) Cf ,U(1) Y q L \ rjK (M 2 • 

°vv -r + <W 73" ) H l \ M v. 



? SU(2) 

>vv 



Cf - (T t 



3 \2 

q~L ' rrXf n/f 2\ , rSU(2)<^A rrN/Ti^ 



c A 



H± (My) + 8vv tttHi (M\ 



2s 



v) 



-X v TlY qL H?(M v ) 



+ 



~2 H l \ M W, 



(73) 
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The unpolarized Born contribution is given in ( 1321) . the counterterms 5C A and 5C N 
are presented in Sect. I3.3[ and the coupling factors are specified by (128"]) and ( 1421) . 
The functions H\(M V ) represent the contributions resulting from the loop diagrams 
of Fig. H They are related to the tensors SA^My) in (SZJ) (J3HD by 



Hl(M^) = -Tr 



, 5 Aii (Mi 



(74) 



The couplings associated with and are the same for q = u and q = 

d. Thus the crossing and CP symmetry relations f[2~Tl) and (1221) imply that these 
functions are symmetric with respect to the transformation t <-> u. In contrast, if^ 
is antisymmetric with respect to t «-> u exchange since the corresponding coupling 
is proportional to T 3 L and has thus opposite signs for q = u and q = d. The 
functions H\(M V ) are presented in Appendix [E] as linear combinations of scalar 
loop integrals. We note that, in contrast to the definition adopted in the case of 
Zj production [9], here we do not include the contributions of the fermionic wave- 
function renormalization constants in Hi (M v ). 

For the IR-singular part of the renormalized one-loop correction we obtain 



qq'—fWg\2 



IR 



^Re 
2tt 



-QqQq'fl + &Qqf2 ~ vQq'fz 



T(l+e)Y}M 



qq'^W°g\2 
I ■ 



(75) 



The IR-singular functions fj R in the MR and DR schemes are presented in Sect. 



4 High-energy limit 

In this section we provide compact analytic expressions that describe the behaviour 
of the IR-finite part of the virtual electroweak corrections in the limit Myy/s — > 
with t/s and u/s constant. In this limit, which is applicable for transverse momenta 
of (9(100 GeV) or beyond, the electroweak corrections are dominated by logarithmic 
contributions of the type ln(s/M^). In Sect. 14.11 we present the asymptotic expan- 
sion of the one-loop corrections, including leading and next-to-leading logarithms, 
as well as terms that are not logarithmically enhanced at high energies. In Sect. 14.21 
we present the two-loop corrections to next-to-leading logarithmic accuracy. 

4.1 Next-to-next-to- leading approximation at one loop 

In this section we discuss the high-energy behaviour of the IR-finite part of the 
one-loop corrections to the qq' — > W a g process, obtained by subtracting the IR 
divergence (j75p from the renormalized one- loop result (ITS"]) . 

j2\ M %7 wa9 \ 2 = El-Mf^'l 2 - El^&^'l 2 - ( 76 ) 
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In the following we present explicit asymptotic expressions for the unrenormalized 
loop contributions, i.e. for the IR-fmite parts ifj' fin of the functions H\ in ( 1731) . 



Using the general results of Ref. [28] , we evaluate the functions H{' &n to next-to- 
next-to-leading logarithmic (NNLL) accuracy. This approximation accounts for all 
contributions that are not suppressed by powers of M^/s. It includes double and 
single logarithms as well as terms that are not logarithmically enhanced in the high- 
energy limit. To simplify non-logarithmic functions of the ratio Mz/My/ we have 
performed an expansion in = 1 — M^/Mf , keeping only terms up to the first 
ordei0 in s 2 ^. The NNLL expansion of if{' fin (My) has the general form 



Hl M {M' } 



NNLL 



V J 



Rc 



i 2 + u 2 
tu 



+ g\{M 2 v . 



i 2 -u 2 

tu 



(77) 



It involves the rational function {i 2 + u 2 )/iu, which has the same angular behaviour 
as the squared Born amplitude (1321) in the high-energy limit, and two other rational 
functions, which describe different angular dependencies. The functions g\ consist of 
logarithms of the kinematical variables and constants. The loop diagrams involving 
Z and W bosons, with mass My = M z , M w , yield 



v) 



9^M 2 ) 



v) 



g?(M 2 ) 

92(M 2 V ) 
9o(M 2 ) 
9?(M 2 V ) 




2+^ UV (M, 



2 

V)i 



-92(M 2 W ), 



0. 
0. 



J In practice we find that all terms of 0{s% l ) cancel in the result. 
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= -2 


21n (4 


gI(M v ) 


= ln 2 ( 




9i(M v ) 


= o, 




gI(M 2 v ) 


= 2 In 


(£)■ 




My J ' 



where A uv is defined in 
we have isolated the term 



(78) 



and, in order to facilitate the comparison with Ref . [H] , 



A,UV 

9o 



(M 



Auv + In 



My, 



(79) 



If we included the fermionic wave-function renormalization constants in the defini- 
tion of the function H±, as we had done for the case of Zj production in Ref. [9] , this 
term would cancel and the function g^My) would be identical to the one obtained 
in Ref. [9]. 

For the loop diagrams involving photons (My = M A ), after subtraction of the 
IR-singular parts, we obtain 



^o N (Mi)=^(M^)-^ + ^= 



go(M 2 A ) 

9l(M 2 A ) 

g\(M 2 A ) 
gl(M 2 A ) 



go(M^) 
gl(M 2 w ) 
g\(M 2 w ) 
g\{M 2 w ) 



7T 



9 

2 

j 

for 
for 
for 



2tt 
73' 



I = X,Y, 

1 = A,N,X,Y, 

I = A,N,X,Y. 



JO) 



The contribution of the counterterms SC A and £C N to the IR-finite part of the 
renormalized result (1TB"]) is obtained by subtracting from (I5U1) the IR-divergent part 
of the wave-function renormalization constants (1621) . This contribution, consisting 
of on-shell self-energies and their derivatives, does not depend on the scattering 
energy. Therefore we evaluate the IR-finite parts of the counterterms in numerical 
form without applying any approximation. Using the input parameters specified in 
Sect. [7] we obtain 



uv 



+ 5.57 x 10" 



5C 



N.fin 




5C 



uv 



- 1.49 x 10" 



The UV divergences 5C A,N |uv [see ( 1601) ] cancel against the Auv-terms in (178 l) -( l80l) . 

The results ( I78|) - (l80|) . for the qq' — > W a g process, are valid for arbitrary values 
of the Mandelstam invariants and can easily be translated to all other processes in 
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(j2J) by means of the relations ( |2~TT) -(l22l). Logarithms with negative arguments in 
( 178l -( l80l) are defined through the usual is prescription, r — > r + is for r = s, t, u. 

In next-to-leading logarithmic (NLL) approximation, i.e. retaining only double 
and single logarithms that grow with energy, the above results assume a particularly 
compact form. In this approximation the counterterms do not contribute, 



sc A N ^ L sc N N ^ L 0, 



(82) 



and for the functions H[' Rn (My), neglecting logarithms of M z /M w , we obtain 



A,fin/ n/T 2 \ NLL 



H?' an (M, 



Vj 



t 2 + u 2 




(83) 



for V = A,Z,W. We note that, owing to H? M (M Z 



NLL 



Hf M {M A ) and 



the NLL contribution of the function H*' n (My) cancels in (1731 . Thus the NLL 
corrections ( l83l) are proportional to the rational function (t 2 +u 2 )/tu, which describes 
the angular dependence of the Born cross section. 



4.2 Next-to-leading logarithms up to two loops 



Let us now present our results for the NLL asymptotic behaviour of the electroweak 
corrections up to two loops. For a discussion of the calculation we refer to Ref . [8] , 
where the same class of corrections has been computed for Zj production. The 
results have been obtained in the M z = M\y approximation. As in the previous 
section, we present results for the IR-finite part of the electroweak corrections, ob- 
tained after subtraction of IR singularities. As discussed in Sect. 13.4} at one loop 
this subtraction is performed in such a way that, to NLL accuracy, the IR-finite 
part corresponds to the complete electroweak correction regularized with a fictitious 
photon mass Ma = M w . The same prescription is adopted at the two-loop level. 

The unpolarized squared matrix element for qq' — > W a g, including NLL terms 
up to the two-loop level, has the general form 



l7r 2 aa s (N 2 - 1) 



V + u 2 
tii 



A (0) 



2tt/ 



+ 



a 
2^ 



34) 



The Born contribution reads 



A {o) 



2s 2 ' 



35) 
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At one loop, the NLL part consists of double- and single-logarithmic terms and reads 



2s 2 



c 



(86) 



Here we used the shorthand L~ = ln fc (|f|/M^ / ) for the logarithms and = 
YqJ (40^) + Cf/s^ are the eigenvalues of the electroweak Casimir operator for left- 
handed quarks. This expression is consistent with the process-independent results 
of Ref. [14J as well as with the NLL part of the one-loop asymptotic expressions 
presented in Sect. 14.11 At two loops we obtain 



NLL 




vf) 



where b\ = —41/(6*^,) and 6 2 = 19/(6^^) are the one- loop /3-function coefficients 
associated with the U(l) and SU(2) couplings, respectively. The LLs as well as the 
angular-dependent subset of the NLLs in (|87|) . i.e. all contributions of the form Lj 
with f = s,t,u, have been derived from Ref. [15]. There, by means of a diagram- 
matic two-loop calculation in the spontaneously broken electroweak theory, it was 
shown that such two-loop terms result from the exponentiation of the correspond- 
ing one-loop corrections. The additional NLLs of the form L| in f)87p have been 
obtained via a fixed-order expansion of the process-independent resummed expres- 
sion proposed in Ref. [IE]. This resummation [16] relies on the assumption that 
effects from spontaneous breaking of the SU(2)xU(l) symmetry can be neglected 
in the high-energy limit. 

Our NNLO predictions include only the LL and NLL terms. Thus they are 
affected by a potentially large theoretical uncertainty, due to missing subleading 
contributions of order a 2 ln fc (s/M^) with k = 2, 1,0. For four-fermion scattering it 
was found that, at s ~ 1 TeV 2 , the two-loop logarithmic expansion has an oscillating 
behaviour characterized by large cancellations between leading and subleading terms 
[3]. In this case the subleading terms play a very important role and the NLL 
approximation yields misleading results. In contrast, in the case of Wj production, 
the relative weight of the LL, NLL and NNLL contributions at one loop indicates 
a fairly good convergence of the logarithmic expansion. Indeed, as can be seen 
from our numerical results in Sect. U\ the one-loop corrections are clearly dominated 
by the negative LL contributions, while the NLL terms are relatively small and 
the NNLL contributions almost negligible. A similar convergence is expected also 
at two-loops, owing to the exponentiation property of the logarithmic corrections. 
Thus our NLL two-loop predictions can be regarded as a plausible estimate of the 
size of the two-loop electroweak effects in high-pp M^-boson production. 
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Figure 4: Tree-level diagrams for the process qq' — > Wgj. 



5 Real corrections 

In order to cancel IR singularities from the virtual-photon corrections, real emission 
corrections need to be calculated. As discussed in Sect. 12.3} all relevant partonic 
reactions are related to the qq' — > W a g^ process through crossing and CP symmetry. 
The tree-level diagrams for this process are shown in Fig. HI 

The squared matrix element for the 2 — > 3 process ab — > W 7 ^, summed over 
polarization and colour as well as averaged over initial-state polarization, can be 
written in a general form 



Ana 



-Q q Q q ,Ht> a \s,t,u,t',u') 



2,ab, 



s, t, u, t' , u') 



aQ q ,H?' ab {s,i,u,i',u') 



with the kinematical invariants defined in Sect. 12.2] In the limit of soft and/or 
collinear photon emission, the squared matrix element (ISSj) exhibits IR singularities. 
To combine these singularities with those originating from virtual corrections we 
have to extract them in analytic form. This is done with the help of the dipole 
subtraction formalism [TTl H8l [T9] . Within this framework the partonic differential 
cross section can be schematically written as 



^ab^Wk-y 

dp T 



Mob / d$ 



M s a u 6 b ($ 3 



+ 



with Nab given in (171). The quantity M reads 



,w " k rFo,sm 



19) 



(90) 



The auxiliary function M^ b is chosen such that it has the same singular behaviour as 
M afe in the soft and collinear limits. This ensures that the difference M ab — can 
be integrated numerically. To compensate for the subtraction, the integral of the 
auxiliary function M^ b , denoted here da^f/dpr, is then added back. The analytical 
form of der^y dpx is obtained after performing integration over the subspace of the 
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radiated photon. The result of this one-particle subspace integration contains sin- 
gular contributions which must be combined with those in virtual corrections. The 
algorithms for constructing the auxiliary subtraction function and its integrated 
counterpart have been developed both for the case of photon radiation off massless 
or massive fermions [T7] and QCD radiation off massless [TH] or massive partons [19J. 
In Sects. 15. II and 15.21 we discuss the application of both formalisms to calculate the 
0(a) real corrections to the Wj production process. In both approaches we use 
expressions for the emission off a massive fermion to describe the emission off a W 
boson, since only soft singularities are present in this case and they depend only on 
the charge of the external particle and not on its spin. 

After adding the real and virtual corrections, collinear singularities remain. 
Final-state singularities are avoided by recombining collinear photon-quark con- 
figurations as discussed in Sect. 12.21 Initial-state singularities are absorbed in the 
definition of PDFs using the MS scheme. 

5.1 Mass regular izat ion 

The formalism of Ref. p2] employs small photon and fermion masses to regularize 
soft and collinear singularities. The subtraction term for the squared matrix element 
is constructed from the appropriate dipole factors. Keeping the original notation of 
Ref. [17J we can write for the process ab — ► W a k^ (where a (b) can be q, q', g) 



{a'=a, b'=b} 



M3b(* 8 ) = -4^ £ I ( Q a Q b 9^r(Pa,P b ,Py) Mf'($ 2 , a 

- QaO gw, T (Pa,PW,Pj) M$ v ($ 2>aW ) 

~ °Qa 9w b aAPW,Pa,P 7 ) Mq' 6 ' ($2,Fa) 
" QaQk g S a tr(Pa,Pk,Pj) M$" (<£ 2 , afe ) 

- Q k Q a g^ T (Pk,Pa,P,) MS' 6 '(<£ 2|to ) + (a<-6) 

+ oQ k g$l >T (pw,Pk,P>r) M a b {®2,wk) j , (91) 
with 

K(®2,nm) = Y,\M^ W ° k (®2,n m )\ 2 Fo,2($2,n m ) • (92) 

Due to Q g = the dipole terms with gluon indices do not contribute to ( I9~TT) and 
for each subprocess the subtraction term M^ b is constructed from six dipole terms, 
characterised by the g snh functions. Expressions for these functions are taken directly 
from Ref. [T7j. In Appendix [F] (see Table [T]) we list all the functions which are used 
to calculate (I9T1) . together with the corresponding equation numbers in Ref. [PT] . 
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For each subprocess the six dipole terms fall into three groups, each containing 
two dipole terms and coming with a specific charge combination, either —Q q Q q > 
or aQq, or —aQ g >. The subtraction term M^ b has then the same structure as M ab 
in (j88p and the IR-singular part of the virtual corrections (1751) . Thus the cancellation 
of singularities can be analyzed for each charge combination separately. 

The construction of the reduced phase space $2,nm follows the prescriptions 
of Ref. [17]. Generally $2,nm is a mapping from the 3-particle phase space into 
a 2-particle phase space. The mapping respects all mass shell conditions. For differ- 
ent types of dipoles, different mappings are necessary. In Table [T] we list numbers of 
equations in Ref. [T7j which we used to perform mapping for the dipole terms appear- 
ing in our calculations. In particular, the observable-defining function Fo }2 (^2,nm) 
in (JS2D is then 

FoA®2,nm) = S( PT - pT,w)0(pT,j ~ p£?) , (93) 



with p TjW and p T;j belonging to $ 2 ,nm- 

The expression for the subtraction term integrated over the phase space of the 
photon reads 



daf 



a I 

"2I ^ a ^ b 



Git (fab] 



a °fwd 



dpi 



+ 



Cdx \g!f(f ab ,x) 

Jo 1 ' 



yS,PT) 



QO fwd 



-Q a cr 



+ 



sub 
1,1 

1 



Gij?W[\?aW 

dx 



r \z.ab^W a k 
acr fwd 



dp T 



dp T 

[s,Pt) 



xs,pt) + (a+->b) 



J ±.ab-*W a k 
UO fwd 



dpi 



xs,p^) + (a*->b) 



-QbV 



GTFM^bW, 



j^ab^Wk 
aa fwd 



dp T 



s,Pt) 



+ 



dx 



acr fwd 



dpi 



xs,pt) + (a<r+b) 



QaQk 



Glf(r ak ] 



j^ab^Wk 



dpT 



{s,Pt) 



+ 



dx 







£r F b (f afc (x),x) 



r \±.ab-+W cr k 
aa fwd 

dpT 



xs,pt) + (a<->6) 



-QbQk 



W F Vbkj 



dpT 



s,Pt) 



+ 



dx 



+QkCr 



G!f(r b k(x) } x) 



A^ab^Wk 



dp T 



QO fwd 



dpT 



xs,pt) + (a<-»6) 



(94) 



25 



The relevant invariants in ( 194)) are defined as r a b = (p a +Pb) 2 = s, f^w = (Pk+Pw) 
s, and 



r a w(x) = (xp a -pw) 2 , 
r a k{x) = {xp a -p k ) 2 , 
r n i = r n i{l). 



(95) 



The terms proportional to d&f^ w " k /dp^ in (G2]) represent the contributions orig- 



inating from the forward hemisphere in the 2-particle phase space [see (!T0]) -(jTT ]) ]. 
The (a<->5) terms are the contributions from the backward hemisphere, and 



a<7 fwd 



dpi 



A^ab^W^k 
uo bkwd 

dp T 



(96) 



Note that the first argument of d(jf^ H/<Tfc /dpi in ( J941 directly indicates the x- 



dependence of the actual values of the t, it invariants defined in ([5]). The plus- 
distributions appearing in eq. (1941 are evaluated according to the prescription 



dx 



g sub (r(x),x) 
i 



da 



dpi 



(xs,p T ) 



dx 



o 



'sub , 



r[x),Xj 



([{T (xs,pr) 9{x-f)- g snh {r{l),x) d ° 



dp 



dpi 



yS,PT) 



(97) 



where f 



Pt + Jpt + 



w 



s guarantees the minimal center-of-mass energy 



to produce the final state. The expressions for the functions G suh and Q suh follow 
directly^ from the results in Ref. [TTj . For the functions G sub they read 



Glf(f) 



G I, 



sub i 
FWH 



Glf(f) 



G 



sub / 
F,FMl 



Re (A IR MR ) + In 2 ( 
+ 2 In 2 



M 2 W , 



3 In 



Re(/ : 
+ ln ^1 
+2Li 2 ( 

Re (/; 

Re(/ ; 



IR 

2,MR 



f 



Mir, 



-In 2 



2 2 „ 

-vr 2 + 4 
3 



M 2 W , 



3M 2 



w 



IR 

MR 



+ ln 2 



IR 

2,MR 



In 2 



2Li 2 



— r 



M 2 W 



.2M 2 , 



-3 In 



2Li 2 
+ 



Mw 



2M^ 



r 



f 



7l 2 1 

H h -. 

6 2 



— r 



1 + In 2 



_ 7T 2 1 

~ 3 + 2' 
M 2 W \ 



(98) 



(99) 
(100) 



5 The function G^^„{f ) has been derived from eq. (4.10) in Ref. [T7] by taking the limit of an 



F,FM 

infinitesimal quark mass 
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4b ' Win f ' -i 



2 \Ml,) 2 \Ml, ) 2 VVr + Afl 



r 




+4L 4^ -4Li4/^) + ^-^ + 3. (101) 



The IR-singular functions //mr cancel against those in the virtual corrections, cf. 
Sect. 13.41 The explicit forms of the functions Q snb are 

gff{r,x) =x(x) + \^\\n(^— ) + 21n(l-x)) + l-x, (102) 



1 — x 



4qed, 



G!tu(r,x) = X (x) + In (^-^) + In ((1 - *)(1 - z 1 (r,x)j)\ 

zi(f,x) — 1 / , , 4Mwi \ 

+ ^rV 3 + Zlif - x) ~ (t - am - J + 1 - {m 



2(1 -x 



+ 1 - x, (104) 



G!f(r, x) = \(.r) + ^— { In ( — ] h, ( I , 

3 



with 



and 



^T-^ (io5) 

= V=7 -21n(l-x)-l}, (106) 

where /Uqed is the factorization scale and m stands for the quark-mass regulator. 
The functions xi x ) are singular. These singularities are related to the collinear 
photon radiation off an initial-state quark and are absorbed in the definition of the 
PDFs, yielding the hadronic cross section finite. The procedure bears complete 
analogy to absorbing collinear QCD singularities into the definition of the PDFs. In 
the MS factorization scheme, the redefinition is achieved by replacing [5] 

fh,q(x,^Q CD )^fh,q(x,^Q CD ,^Q ED ) — —Qgf fh,q(~ ■> ^QCD) A*QEd) [x( 2 )]+ ■ 

Z7l Jx Z Z 

(107) 



5.2 Dimensional regularization 

In an independent calculation we used the results of Refs. [TBI [E] to evaluate the 
dipole subtraction terms and their integrated counterparts. The formalism of [T8|fT9] 
is concerned with QCD radiation and expressions for dipoles are given as matrices in 
colour and helicity space. Since we consider photon emission off a fermion line, the 
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colour and helicity structure disappears and the dipole matrices reduce to simple 
expressions. More precisely, to adapt the formalism Refs. [TBI [T9] for the calculation 
of QED corrections, we make use of expressions describing gluon radiation off a 
fermion line in Refs. (TgJ [TU] and replace 



as — ► a, 



o'iQi, C-p — > Ql, Tr — ► 1, Ca 







(106 



where Tj indicates the colour of the emitting parton, Qi is the electric charge in units 
of the positron charge for this parton, and <7j = +1 (—1) for incoming (outgoing) 
partons. Adopting notation analogous to Refs. [TBI IE], the subtraction term for the 
process ab — *> W a k^ can be then written 



M ab 



sub 



($3 



af,b j_ -t-)07 j_ T> a 

QED U W, QED "T u -y 



-yW, QED 



+ u k, QED 



+ ^fc, QED + («^) 



+ D~/k,W, QED + DjW,k, QED ) 



where 



■^F, QED 



F 0t2 (^2,nm) Vf(PW, Pk, P^ Pa, Pb) 



replacements of eq. [|i08|l 



(109) 



110) 



It is understood in eq. (11091) that dipole subtraction terms with a gluon index do 
not contribute to M^ b . In a complete analogy to eq. (JHJJ), for any initial state ab 
the expression for M^ b is constructed from six dipole subtraction terms I^qed) 
each associated with one of the three possible charge combinations —Q q Q q >, crQ q or 
—aQ q i. The dipole subtraction functions T> F are taken directly from Refs. [TH1 [TH] . A 
list of the functions Dp used to calculate the subtraction term M^ b in ( 1109[) . together 
with the corresponding equation numbers in Refs. [THl UH], is presented in Table El 
Appendix [Fl Additionally, for each dipole subtraction term appearing in ( 1 109ft we 
include a description of its type. The mappings from $ 3 to $ 2 ,nm agree between the 
formalism of Refs. [TBI US] and [T7j. However, for the sake of completeness, Table El 
contains numbers of equations which provide mapping formulae in Refs. [TBI [T9] . 
The function Fo,2 in (HI 01) is given by expression (1931) . 

Moreover, apart from the final-state emitter, final-state spectator case, i.e. the 
dipoles T> k W q ED and T> Wk QED , there is a direct correspondence between the 
dipole subtraction terms in the two formalisms of the form 



'„ab) 



iW , 



V$SD -QaQb 4™ £ g^ T ( Pa ,P b ,P,)K(ha 

T=± 

^ QED QaCT 47m J2 9atA^PW,P,)K^2,aX 

T=± 

^V,QED Qa^ 47m Yl gWaAPWiPa^^Kihwa), 

T=± 

VPqfSD QaQk 4^ £ gtr(Pa,Pk,Pj)K b ^2,ak), 

T=± 

^ fc ,QED QaQk 47TC £ f/g T (pk,Pa,Ih)K(&V 



:,ka) 



'111) 
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The subtraction term integrated over the photon phase space is constructed 
according to 



daf _ a 
dpT 2rr 



I ab (S,PT) + £dx (K, ab (x,s,p T )+V ab {x,s, PT/ 



(112) 



where the expressions for X, /C and V functions follow from results for the integrated 
dipole functions in Refs. [18, 19J after performing replacements of eq. fllOSp . For the 



photonic corrections to any of the subprocesses ab — > W a k we can write 



l ab (s,p T ) =Q a Q b 

- QaCT 

- Qb<? 
QaQk 

- QbQk 

+ Qk<y 



1(Sr, 



A^ab^Wk 

LIU r 



fwd 



1'{s a 



dp-] 
, da 



s,Pt) + (aW>) 



w , 



ab^Wk 
fwd 



dpi 



s,Pt) + 



Til* \ uo fwd 
X (S bW ) — 

dp T 

A-ab^Wk 
i{Sak)^f (S,PT) + 

dp T 

i(hk)^ — (s,pt) + 



dpi 



T (skw 



aa fwd 



dpi 



s,Pr) + 



(a<->6) 
(a<->-6) 
(a<->6) 
(a<r->b) 
(a<->6) 
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where s nm = 2p n p m , and we make use of Q g = 0. As in (I'M)) , the terms proportional 
to d<jf^ M/<Tfe /dpT originate from the forward hemisphere and the (a<-+b) terms from 
the backward hemisphere. The integrated dipole functions in (I113p read 



2\ 8 



4vr 2 



r(l + £ )Re(A IR DR ) 



-In 



2 / °nm 



M W, 



+ 3 In 



M w, 



2\ £ 



In 



r(l + e)Re(/ : 
In 



IR 

2,DR 



-In 



2 / °nm 



1 urn 



+ ln 



M 2 W 



Myy + <S n?Tt/ 



-I In 



In 



14^ ^ °nm / 
Sum 



+ 



) + 2Li 2 



+ 3 In 1 



\M W + 



3M, 



7T 2 - 6. 



(114) 
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The structure of the singular terms //dr. i n pi4p is kept the same as in ( 175|) to 
manifestly show cancellation of singularities between virtual and real corrections. 
For the x-dependent functions we have 



}C ab (x,s, PT )= (Ql + Ql 

+ 2QaQb 
- Qa(? 



K{x 



A~ab^W a k 
QO fwd 



dpi 



j -ab^W^k 

jC(x)=^f (xs,pr) + {a^b) 



{ i^.ab^W a k 
w/ - \ ua fwd 
/C [X, S aW ) 1 

dp T 

^frab^W^k 

IC'(x, s b w) — "ir (xs,pt) + (a<r+b) 



xs,pt) + (a<->6) 



— (QaQk + QbQk) 



dpT 

A-ab^W rT k 

K!'(x) fwd 



xs,p-r) + (a<->6) 



V a \x } s,pr) 



ZQaQb 

- QbO 

QaQk 

— QbQk 



V(x, s ab 



Q(J fwd 



dp 



'T 
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^ab-*W a k 
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dpT 
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Note that in contrast to eq. (II 13j) . the quantity s nm in (j!15p can be implicitly 
dependent on the fraction x. More precisely, it is the case if s nm involves the 
momentum of a final-state particle. The final-state momentum belongs then to the 
phase space for which the squared center-of-mass energy is xs = 2xp a pb [19|. 

The evaluation of the terms involving the plus-distribution is carried out as 
indicated in Ref. [19], i.e. according to 



o 



dx [TZ(x, s nm (x))} + (xs. i> }: ) 



dpi 



dx 



da 



ft(x } s nm (x))- — (xs,p T ) 6(x - f) 
dp T 



fc{x, s nm (l)) 



da 
dp T 



[s,Pt) 



'1181 



with f = ( p T + ^p\ + M^j I s. 

In the formalism of Refs. [HI [19] the collinear counterterms associated with PDF 
renormalization are included in the expressions for integrated dipole functions, i.e. 
the final results which we use are free from collinear singularities. The expressions 
presented here are calculated using the MS factorization scheme. 

As can be seen from the presented formulae, the explicit expressions for the 
integrated dipole functions in the two formalisms are different. In particular, the 
expressions for the end-point contributions have different forms due to specific con- 
ventions wrt. calculating the plus-distribution terms in the two formalisms. However 
we have checked that, after subtraction of the IR singularities, for each charge com- 
bination apart from aQk the integrated dipole contributions to dcr^/dpr i n t ne two 
formalisms are equivalent. 



31 



6 Checks 



Every part of the presented calculation has been performed in two completely in- 
dependent ways. The algebraic reductions were done using two different Mathe- 
matica [30] codes. For the numerical evaluation we have implemented the results 
in two independent Fortran programs. Comparing the results at numerical level 
we find agreement within the statistical errors. 

Furthermore, in order to control the correctness of our results we performed 
various consistency checks. On the side of virtual corrections we have verified that 
the one-loop corrections (1481) satisfy the Ward Identity 

eVMVguvfa) [5A^(M*)u x ] u(p q ) = for I=A,N,X,Y. (119) 

A similar Ward identity holds for the lowest-order amplitude^). The cancellation of 
the ultraviolet divergencies has been verified analytically and numerically. For the 
numerical evaluation of the loop integrals we use a set of routines by A. Denner 
and, alternatively, the FF library [31]. The NLL approximation that was derived 
from the full one-loop calculation, has been checked against results from the general 
derivation of NLL terms [II] . Also the IR-singular contributions in the high-energy 
limit have been reproduced within this framework. 

The squared matrix element for the real corrections was checked numerically 
against MadGraph [32]. The cancellation of IR singularities between real and vir- 
tual corrections was done analytically using the dipole formalism. The subtraction 
terms were derived and implemented in two different ways, using the mass regu- 
larization of IR singularities and the dimensional regularization. The phase-space 
integration for the real corrections was performed with adaptive Monte-Carlo inte- 
gration using VEGAS [33]. Detailed comparisons at analytical and numerical level 
were performed, and the agreement between the predictions generated within two 
different regularization schemes provided a strong check on the calculation of the 
real corrections. 



7 Numerical results 



In this section we present numerical predictions for the large-pr production of W 
bosons at the hadron colliders LHC and Tevatron. The following input parame- 



6 We note that the abelian one-loop contribution satisfies two additional Ward identities 
PWii£* v (Pg)v(Pq) 5A^ A (My)uj\ u(p q )=0, pw^PgvViPq) ^i,a( M v-) w a u(p g )=0. 

Similar identities for the N-, X- and Y- form factors exist but are less trivial due to the non- 
vanishing contributions from would-be Goldstone bosons on the right-hand side. This means that 
the calculation of the unpolarized cross section requires the use of the exact expression for the W- 
boson polarization sum. Instead, owing to (|119[) . the gluon polarization sum can be implemented 

as ~9vv' ■ 
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ters are used: = 1.16637 x lO^GeV" 2 , M w = 80.39 GeV, M z = 91.19 GeV, 
m t = 171.4 GeV, Mh = 120 GeV. For the numerical values of elements in the CKM 
quark mixing matrix we refer to [34] . 

The hadronic cross sections are obtained using LO MRST2001 PDFs (35] • We 
choose /iq CD = p\ as the factorization scale and, similarly, as the scale at which 
the strong coupling constant is evaluated^- We also adopt, in agreement with the 
value used in the PDF analysis, the value as(Mf) = 0.13 and use the one- loop 
running expression for tts^Qco)- I* 1 our calculations of the real corrections we 
choose the MS factorization scheme with the scale /iq ED = M^- We note that in 
order to consistently include 0(a) corrections in a calculation of a hadronic cross 
section, PDFs that are used in the calculation need to take into account QED 
effects. Such PDF analysis has been performed in [36J and the 0(a) effects are 
known to be small for /iQED ^ 100 GeV, both concerning the change in the quark 
distribution functions (below (9(1%) [37]) and the size of the photon distribution 
function. Moreover, the currently available PDFs incorporating 0(a) corrections, 
MRST2004QED [36], include QCD effects at the NLO in as- Since our calculations 
are of the lowest order in QCD, and QED effects on PDFs are estimated to be small 
for /iQED 100 GeV, we prefer to use a LO QCD PDF set without QED corrections 
incorporated, rather than MRST2004QED, and we set /iqed = MyJo- Moreover we 
do not include photon-induced contributions, which are parametrically suppressed 
by a factor a /as- However, in the concurrent to this paper (and subsequent to 
Ref. pi]), work of Ref. [22], it has been reported that photon- induced contributions 
are of numerical significance for large pr W^-boson production at the LHC. Estimates 
of the exact size of these effects are obscured by large theoretical uncertainty on the 
photon's PDF, as demonstrated in Ref. |22j . 

We choose the following values of the pp-cuts: p™™ = 100 GeV for LHC and 
p™j = 50 GeV for Tevatron. The value of the separation parameter below which 
the recombination procedure is applied is taken to be R sep = 0.4. The dependence of 
our predictions on R sep is negligible. We have verified that the shift of the transverse- 

7 Note that when calculating the contribution to the hadronic cross section coming from the 
subtraction term in the real corrections, we take the transverse momentum of the W boson in the 
reduced phase space, Pt,w, as the factorization scale and the argument of as- 

8 The use of different factorization scales, /iQCD = Pt and /iQED = Mw , is due to the fact that 
Mqcd and /iQED play a different role in our calculation. The dependence on /iQCD is due to the LO 
evolution of the PDFs and represents an effect of 0(as ln(/iQCD / Mo)), where /io is the scale at which 
the PDF evolution starts. This dependence would be compensated by NLO QCD contributions of 
0{as ln(pr/MQCD)) and, although QCD corrections are not included in our calculation, choosing 
MQCD = Pt we can absorb large NLO QCD logarithms of the scale pr in the LO PDF evolution. In 
contrast, the /iQED dependence of our predictions is due to 0(a ln(pi>//iQED)) terms in the photon 
bremsstrahlung corrections. This dependence is not compensated by the PDF evolution since we 
use a PDF set that does not include QED effects, assuming that these effects are negligible. This 
approach makes sense only if the scale /iQED is chosen in such a way that the (potential) impact 
of QED effects on the PDFs is very small. In Ref. [37J it was shown that the QED corrections to 
the PDFs grow with /iQED but do not exceed one percent for /iQED ^5 100 GeV. This motivates 
our choice /iQED = M\y for the QED factorization scale. 
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momentum distribution induced by variations of this parameter in the range 0.1 < 
-Rsep < 1-0 does not exceed a few permille. 

Our lowest-order (LO) predictions result from (1321) . The next-to- leading order 
predictions (NLO) include the LO+virtual contributions (173"]) and the real brems- 
strahlung (1591) . We also study the relative importance of the IR-finite parts of the 
virtual (NLO V i rt ) and real (NLO rea i) contributions. These IR-finite parts are con- 
structed by subtracting the IR divergence (175!) from the virtual corrections and 
adding it to the real ones. The next-to-leading-logarithmic (NLL) and next-to- 
next-leading-logarithmic (NNLL) predictions^ are obtained adding to the LO the 
approximations (1HB"]) and (1771) for the NLO virt part and neglecting the NLO rea i part 
of the corrections. As we will demonstrate, for the case of fully inclusive photon 
radiation neglecting this piece provides a good approximation of the complete calcu- 
lation. The next-to-next-to-leading order predictions (NNLO) include the full NLO 
results plus the two-loop NLL corrections (1571) . 

The LO transverse-momentum distributions for pp — > W + j and pp — > W~j at 
the LHC are shown in Fig. EK- In Fig. [5b and Fig. [5fc we plot the relative size of the 
NLO, one-loop NLL, one-loop NNLL and NNLO corrections wrt. the LO predictions 
for W + and W~ production, respectively. The behaviour of the relative corrections 
to W + and W~ production is very similar. As expected, the importance of the 
NLO contribution increases significantly with px and leads to a negative correction 
ranging from —15% at px = 500 GeV to —43% at p T = 2 TeV. We also observe that 
the one-loop NLL and NNLL approximations are in good agreement (at the 1-2% 
level) with the full NLO result for px > 100 GeV. The difference between NLO and 
NNLO curves is significant. The two-loop terms are positive and amount to +3% 
at pt = 1 TeV and +9% at px — 2 TeV. This shifts the relative corrections for W + 
production up to —25% at px = 1 TeV and —34% at px = 2 TeV. 

The IR-finite parts of the virtual (NLO v i rt ) and real (NLO rea i) corrections to 
W + production at the LHC are shown separately in Fig. Eh- Fig. [6b shows the 
relative size of the NLO V i rt and NLO rea i corrections wrt. the LO predictions. The 
NLO v i rt contribution dominates the full NLO correction and amounts up to —42% 
at px = 2 TeV. The NLO rea i part contributes with a smaller and nearly constant 
correction of about —1% in the entire px-range. This means that, for the case of 
fully inclusive photon radiation, the NLO vrr t part represents a good approximation 
of the full NLO correction. 

The high-energy behaviour of the NLO v i rt part is described by the compact NLL 
and NNLL approximations presented in Sect. HI The quality of these approximations 
is shown in Fig. [7J We observe that the NLL approximation works well differing from 
the exact NLO v i r t result by less than 1% for px > 200 GeV. The quality of the NNLL 
approximation is of the order of one permille or better in the entire px-range. 

For less inclusive observables where a veto on hard photons is imposed, the 
NLO rca i contribution can become important. Fig. [H] shows the relative NLO rea i cor- 

9 For details concerning the treatment of angular-dependent logarithms at the NLL level we 
refer to Ref. 0. 



34 



10 1 




LO 


10° 






> 10" 1 

0) 






LD in? 

\ 10 






Q. 1f)-3 






10" 5 






10" 6 


- (a) 




0.00 


a] 


NNLO/LO 1 


-0.10 


1 ■ 1 


NLO/LO 1 
NLL/I 1 

1 M l_ l_ / l_ V_/ ± 


-0.20 


' ' 1 " f _ 1 


NNLL/LO - 1 ■ 


-0.30 






-0.40 


(b) W+ 




0.00 




NNLO/LO 1 


-0.10 




NLO/LO 1 
NLL/LO 1 


-0.20 




NNLL/LO -1 ■ 


-0.30 






-0.40 


- (c) W" 





200 400 600 800 1000 1200 1400 1600 1800 2000 

Pr [GeV] 



Figure 5: Transverse-momentum distribution for V^-boson production at the LHC. 
(a) LO distribution for pp— >W + j and pp— >W~j. (b) Relative NLO (dotted), NLL 
(thin solid), NNLL (squares) and NNLO (thick solid) electroweak correction wrt. the 
LO distribution for pp^W + j. (c) Relative NLO (dotted), NLL (thin solid), NNLL 
(squares) and NNLO (thick solid) electroweak correction wrt. the LO distribution 
for pp^W~j. 
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Figure 6: IR-finite parts of the virtual (NLO virt ) and real (NLO rea i) contributions 
to the pT-distribution of W bosons in the process pp^W + j at ^/s = 14 TeV. 
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Figure 7: Relative precision of the high-energy approximations at one loop in the 
process pp^W + j at yfs = 14TeV as a function of pt'- NNLL (solid) and NLL 
(dashed) wrt. the IR-finite part of the exact one-loop result (LO+NLO virt ). 
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Figure 8: Relative size of the real correction (NLO rea i) wrt. the LO for fully inclu- 
sive photon radiation and the case where visible photons are rejected, plotted as a 
function of pp for the process pp^W + j at y/s = 14 TeV. 
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rections for W + production. We compare the fully inclusive photon radiation with 
the case where visible photons with px, 7 > lOGeV and R(j,j) > 0.4 are rejected. 
This veto leads to a significant enhancement of the (absolute size of the) NLO rea i 
part, which can exceed —5% for px > 1 TeV. 

To underline the relevance of the large electroweak corrections for W produc- 
tion at the LHC, in Fig. [9^ and Fig. [9b we present the relative NLO and NNLO 
corrections to the W + and W~ cross sections integrated over px starting from 
Pt = Pt u *5 as a function of p T ut . This is compared with the statistical error, de- 
fined as A(j stat /(j = l/y/N with N = £ x BR x o~ lo . The branching ratio BR = 2/9 
accounts for the full efficiency of jy-detection in the ez/ e and fii/^ modes (for this 
estimate we ignore experimental efficiencies and cuts) and we assume a total inte- 
grated luminosity £ = 300fb _1 for the LHC [38]. It is clear that the size of the NLO 
corrections is much bigger than the statistical error. Indeed, already for £ = 3fb _1 
and pt < 800 GeV they correspond to a two standard deviation effect. Also the 
difference between the NNLO and NLO corrections, due to two-loop logarithmic 
effects, is significant. In terms of the estimated statistical error, these two- loop 
contributions amount to 1-3 standard deviations for px of (9(1 TeV). 

Ratios of px-distributions for W + , W~ , Z bosons [9] and photons [ID] , in con- 
trast to the distributions themselves are expected to be relatively insensitive to QCD 
corrections and theoretical uncertainties associated with as and PDFs. These ra- 
tios lead to important experimental tests of W and Z couplings in the high-energy 
region. For W + and W~ the ratio is presented in Fig. [TOa . The LO value increases 
from 1.5 at px — 100 GeV to 3.4 at px = 2 TeV. As already observed, the (rela- 
tive) electroweak corrections to the W + - and VT~-boson px-distributions are almost 
identical. In consequence, the LO, NLO and NNLO curves in Fig. ITDa overlap. In 
contrast, the impact of the electroweak corrections on the W + ratio (Fig. [TDb) at 
the LHC is clearly visible. The LO prediction, ranging from 1.4 to 2.5, receives a 
negative NLO correction that grows with px and amounts to —0.5 for px = 1 TeV. 
At px = 2 TeV the difference between the NNLO and NLO curves is about 0.2. 

The ratios of px-distributions for W + /Z and W~/Z are shown in Fig. [TTa and 
Fig. [TTb . respectively. For the W + jZ ratio the LO prediction ranges from 1.5 to 2. 
For px > 1 TeV it is reduced by 0.09 to 0.18 by the NLO electroweak corrections. 
The logarithmic two-loop corrections are small. A qualitatively similar behaviour is 
observed for the W~ jZ ratio. 

The results of a similar analysis for W + production at the Tevatron (-^/s = 2 TeV) 
are shown in Figs. [T3?HT5l (the px-distributions for W + and W~ production are ob- 
viously identical). The LO px-distribution is shown in Fig. fT2"k. the relative NLO, 
NLL, NNLL and NNLO corrections in Fig. fT2b . The NLO corrections grow with 
px and reach —11% at px = 400 GeV. The one-loop NNLL and NLL approxima- 
tions describe the exact NLO results with 1% and 3% precision, respectively. The 
dominant two-loop effects have little impact on the size of the corrections. 

The quality of the high-energy approximations wrt. the IR-finite part of the 
virtual corrections (NLO v i r t) is shown in Fig. [131 In the px-range under considera- 
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Figure 9: Relative NLO (dotted) and NNLO (solid) electroweak corrections wrt. the 
LO and statistical error (shaded area) for the unpolarized integrated cross section 
for (a) pp—>W + j at y/s = 14 TeV and (b) pp^W~j at \/s = 14 TeV as a function 
of p c ^\W). 
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Figure 10: Ratio of the pp-distributions for the processes (a) pp^W + j and pp—>W~j 
and (b) pp^>W + j and pp—'-'jj at yfs = 14TeV: LO (thin solid), NLO(dotted) and 
NNLO (thick solid) predictions. 
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Figure 11: Ratio of the pp-distributions for the processes (a) pp^>W + j and pp^-Zj 
and (b) pp^>W~j and pp^-Zj at y/s = 14TeV: LO (thin solid), NLO(dotted) and 
NNLO (thick solid) predictions. 
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tion both approximations are less precise than at the LHC, nevertheless the NNLL 
approximation in sufficient for all practical purposes. 

In Fig. [TUthe relative NLO and NNLO corrections to the integrated cross section 
for j»t > Px ut are compared with the expected statistical error for an integrated 
luminosity C = 7fb _1 [39J. The size of the NLO electroweak corrections is above 
the statistical error for a significant range of px-values. Therefore they should be 
included in the analysis when considering precision measurements. In contrast, the 
impact of the dominant two-loop corrections is negligible. 

The effect of the electroweak corrections on the ratios of pp-distributions for 
W/Z and for W/j is shown in Fig. [T5"a and Fig. [Tot ), respectively. 

8 Summary 

In this work the electroweak corrections to large transverse momentum production 
of W bosons at the hadron colliders Tevatron and LHC were evaluated. The con- 
tributions from real and virtual photons cannot be separated in a gauge-invariant 
manner from purely weak corrections and were thus included in our analysis. Soft 
and collinear singularities were regulated by introducing a small quark mass and 
a small photon mass and, alternatively, by using dimensional regularization. The 
real photon radiation was evaluated using the dipole subtraction formalism. The 
agreement between the results derived in the two regularization schemes has been an 
important cross check of the calculation. Numerous additional tests were performed 
to ensure the correctness of the result. 

At the Tevatron, p^- values up to around 300 GeV can be reached with reasonable 
event rates. In this region the 0(a) electroweak corrections reach up to —10% and 
are thus of relevance for precision measurements. Two-loop electroweak corrections 
are negligible at the Tevatron. With below 400 GeV the relative rates for W, Z 
and 7 production are hardly affected by electroweak corrections. 

In contrast, for transverse momenta in the TeV region accessible at the LHC, 
electroweak corrections play an important role. The 0(a) corrections lead to a 
reduction of the cross section by about —15% at transverse momenta of 500 GeV 
and reach more than —40% at 2 TeV. The logarithmically dominant terms were 
extracted from the exact expression of the virtual corrections and agreement was 
found with the predictions based on the process-independent analysis of electroweak 
Sudakov logarithms. If no cuts on real photons are applied, the contribution of the 
real photon emission is numerically small (about 1%) and almost independent of pp- 
Numerically the NLL and NNLL approximations give a good description of the full 
0(a) result with an accuracy of about 1-2%. Considering the large event rate at 
the LHC, leading to a fairly good statistical precision even at transverse momenta 
up to 2 TeV, we evaluated also the dominant (NLL) two-loop terms. In the high-pr 
region, these two-loop logarithmic effects increase the cross section by 5-10% and 
thus become of importance in precision studies. We also studied the relative rates 



42 



for W + , W~, Z and 7 production, which are expected to be stable with respect to 
QCD effects. The electroweak corrections cancel almost completely in the W + /W~ 
ratio. In contrast, their impact on the W + /Z and the W + /7 ratios is significant 
and leads to a shift of 0(10%) for p T > 1 TeV. 
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Figure 12: Transverse-momentum distribution for W^-boson production at the Teva- 
tron. (a) LO distribution for pp^W +{ - ] j (solid), (b) Relative NLO (dotted), NLL 
(thin solid), NNLL (squares) and NNLO (thick solid) electroweak correction wrt. 
the LO distribution for pp— >W + (~> j. 
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Figure 13: Relative precision of the high-energy approximations at one loop in the 
process pp— >W +<K ~^j at y/s = 2TeV as a function of p^: NNLL (solid) and NLL 
(dashed) wrt. the IR-finite part of the exact one-loop result (LO+NLO virt ). 
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Figure 14: Relative NLO (dotted) and NNLO (solid) electroweak corrections wrt. 
the LO and statistical error (shaded area) for the unpolarized integrated cross section 
for pp—*W + ('j at yfs = 2 TeV as a function of p^ t (W). 
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Figure 15: Ratio of the px-distributions (a) for the processes pp-^W +<y >j and 
pp^Zj and (b) for the processes pp— >W + (~> j and pp^'jj at -/s = 2 TeV: LO 
(thin solid), NLO (dotted) and NNLO (thick solid) predictions. 
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A Recombination and exclusive Wj cross section 



As discussed in Sect. 12.2} the recombination prescription that we use to regularize 
photon-quark final-state collinear singularities implies a different treatment of final- 
state quarks and gluons. While for final-state gluons we apply a cut on £>t, 9 within 
the entire phase space, for final-state quarks the recombination effectively removes 
the cut on pr, 9 inside the collinear cone R(q,j) < R scp . As a consequence the 
recombined gq' — > W a q r ) cross section (fI51) has a logarithmic dependence on the 
cut-off parameter -R sep . In order to quantify this R sep - dependence, let us consider 
the contribution of real photon radiation inside the recombination cone. To this 
end, assuming that the cone is sufficiently small (-R S ep "C 1), we adopt a collinear 
approximation 

f da 9q '^ w ^ = a 9q '^ w ° q f 1 dzFgJz), (120) 

J R(q,-y)<R scp JO 

where^l z = PT,-y/(PT, q + Pt,~j) = 1 — PT, q /PT,w is the photon momentum fraction 
and [IE] 

a ®l d ^ (4vr/i 2 ) e /- fc i. m ax dk] 



2tt "' v "' T(l -e)Jo (kj_) 1+E 
' iQ %W#T+^(^ 2 ) (121) 



2tt 97V J eT(l-e) 



with 



aQ 



2 

2tt 



P„Jz) In 



H 2 



A- 2 

l.max / 



(122) 



Here P qi (z, e) = P qi {z) — ez with P qi {z) = [1 + (1 — z) 2 ]/z is the q — > 7 splitting 
function in 4— 2e dimensions, fc^ is the photon tranverse momentum wrt. the photon- 
quark system, and /c_i_ jmax = z(l — z)R sep p TjW . The 1/e collinear singularity resulting 
from inclusive photon radiation, i.e. integrating over the complete energy spectrum 
< z < 1, cancels against the virtual corrections. 

The -R S ep- dependence of the recombined cross section (TTB1) is due to the fact that, 
inside the recombination cone quarks with p^ q < p^ m - (or equivalently photons with 
z > 1 —Ptj/pt,w) are n °t rejected. Thus the variation of <r rec . induced by a rescaling 
-Rsep -> £sep-Rse P amounts to 

a £Uw q = ln ^e P ^ («) with ^ min = 1 - p£*/p T , W . (123) 

For relatively small transverse momenta {pt,w — ^P^ lI j) a rescaling of R sep by a factor 
£sep — 10 shifts the gq' — > W^qij) cross section by less than 2 (0.5) permille for up- 



3 Here we assume lowest-order kinematics, i.e. px,g + Pi,^ — Pt,w m the collinear region. 
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(down-) type quarks. Moreover it is obvious that at highpr.W) where z min — > 1, this 
effect tends to disappear. 

Let us now compare the recombination procedure with a realistic definition of 
exclusive pp — > Wj production, where final-state quarks (a = q) and gluons (a = g) 
are subject to the same cut p^ >a > p^j within the entire phase space (including 
collinear quark-photon configurations). Since the recombination procedure does 
not affect final-state gluons, only channels involving final-state quarks need to be 
considered. The difference between the recombined gq' — > W u q r y cross section ffl"8l) 
and the exclusive cross section (TE?j) corresponds to the contribution of hard collinear 
photons with R(q, 7) < R sep and z min < z < 1. This collinear hard-photon radiation 
can be described by means of quark fragmentation functions [IHJ EH EHJ 06j H3, HHj 
as 

Aa cxcl . = a^ w ^ - d«£? w '*> = a^ w °« f dzV q ,(z). (124) 

" ■^•min 

Here the effective quark fragmentation function T> qi {z) = F q7 (z) + D qi (z) consists 
of the perturbative contribution F qi and the bare fragmentation function D qi . The 
collinear singularity resulting from the perturbative contribution is factorized into 
the bare fragmentation function at the scale fi, such that in the MS scheme [1H] 

V qi (z) = F qi {z, fi 2 ) + D„{z, /i 2 ), (125) 

and the renormalized fragmentation function D qi can be extracted from experimen- 
tal measurements. Using the parametrization [441 148] 

D g7 (z, /i 2 ) = ^ \-P qi {z) ln(l - zf - 13.26] , (126) 

Z7T L J 

obtained by the ALEPH collaboration at /io = 0.14 GeV, we arrive at 

aQ 2 g 



P„(z)i a ( zR ™° PT ' w ) 2 + z- 13.26 



V 



/'0 



(127) 



With this expression we derive a conservative upper bound for A<T exc i.. To this 
end we consider Q q = 2/3, R sep ~ 1, and a wide range of transverse momenta, 
2px'j — Pt,w < 2 TeV. With these parameters we obtain 

< 2 x 1(T 3 . (128) 

a 

We conclude that, for R sep < 0(1), the recombined cross section has a negligi- 
ble dependence on the recombination parameter R sep and provides a fairly precise 
description of exclusive pp — > Wj production at high transverse momentum. 
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B Standard matrix elements 



The algebraic expressions involving external momenta, Dirac matrices, spinors and 
gauge-boson polarization vectors have been reduced to a set of 10 standard matrix 
elements 

Si = v(pq)SruJLu(p q >) e*^p w )el(p g ), (129) 

with 





= 7 ypw 


-fan > 


st 


= {jfw-: 




st 


= Yp u w, 




st 


= -7X> 




^5 


= 7%> 






= -Yp$, 




st 


= (jfar-: 




^8 


= (jfo-: 


i>g)PqP V q >, 




= (jfar-: 




°10 


= (jhr-; 


$g)PqPw- 



These algebraic expressions correspond to the massless subset of the standard matrix 
elements of Ref . 125 1. 



C Scalar loop integrals 

In this appendix we list the scalar loop integrals Jj(M v ) that contribute to fj48l) . 
The symbols Jj are chosen in analogy with Ref. p]. For convenience, to denote 
constant terms we define 

J (My) = 1. (131) 

For the scalar integrals A , B , Cq and D we adopt the notation of FeynCalc [40J. 
However, we choose their normalization according to Ref. [2S]> i- e - we include the 
factor (27Tyu) 4 ~ D which is omitted in the conventions of FeynCalc. 
The UV-divergent one- and two-point functions are denoted as 



Jla{My) 


= B (m 2 ;M 2 ,m 2 ), 


Jlb(My) 


= B (m 2 ;M^,m 2 ) = Jia{Ml 


J 2 (My) 


= B (p 2 w ;m 2 ,m 2 ), 


J3(My) 


= B Q (p 2 w ;M^,M v ), 


Jl{My) 


= B (s;m 2 ,m 2 ), 


MMy) 


= B (u;M 2 ,m 2 ), 
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MM 2 ) 



B Q {u\M\ 



m 



Bo(i; My, m 2 
J 6b (M v ) = B (t;M^,m 2 



JUMw) 



(132) 



The remaining loop integrals are free from UV singularities. The following three- 
point functions are finite if My and the W-boson transverse momentum are non- 
vanishing: 



J 7 {My) 


= C (s,m 2 , 


m 2 ; 


m 2 ,m 2 ,M 2 ), 








UMy) 


= C (u,p 2 w 


m 2 


M v ,m 2 ,m 2 ), 








J9a(My) 


= C (u,p 2 w 


m 2 


m 2 ,M^,M 2 ), 








MMy) 


= Cq(u,p 2 w 


m 2 


m 2 ,M 2 ,M^) -- 


= J9a(My) 




MK) 


— Co(^)PvF) 


m 2 ; 


M v ,m 2 ,m 2 ) = 


UM 2 ) 


t*- 


>u 


Jlla(My) 


— Co^jply, 


m 2 ; 


m 2 ,M^M 2 ) = 


~- JUMy) 




Jllb(My) 


= CoitjPw, 


m 2 ; 


m 2 ,M 2 ,M^) = 


-- MM 2 ) 


i*->u ' 



(133) 



In addition, the box diagrams bl-b3 in Fig. [21 provide the following combinations 
of three- and four-point functions 



MMy) 



D (m 2 , OjPiyi 71,2 j u, s; M v , m 2 , m 2 , m 2 ) 



1 



x 



su + {t 

u-pw)CQ{u,p 2 v ,m 2 ]My,m 2 ,m 2 ) + uC (u,0,m 2 ; M v ,m 2 ,m 2 ) 
+ {s -p H /)C (s,p 2 y ,0;m 2 ,m 2 ,m 2 ) 



{l)My 



Jl3(My) 

Jl4a(My) 



Jub(M 2 ) 



MM 2 ) 



t<-*u 

Doipiy.nr.O.ur.l.a-.M 2 .}! 2 2 



m 2 , m" 



iC {i, 0, m 2 ; My, m 2 , m 2 ) + uC {u, 0, m 2 ; M^r, m 2 , m 2 



tu - tM^y - uMl 



1 v 



M^M V 



JlAa{Ml 



(134) 



For non-vanishing My and PF-boson transverse momentum, the functions Ju^Jub 
are finite. The fact that the scalar four-point functions in f)134p appear always 
in combination with three-point functions is due to the cancellation of the collinear 
singularities that are associated with the gqq vertex [9] . Although these singularities 
are present in individual -Do and Co functions, they always cancel in the complete 
result for box diagrams. 
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D Infrared singularities 



The scalar integrals J{(M V ) in Appendix O contain soft and collinear singularities 
that appear when My = Ma —> and m — ► 0. As discussed in Sect. 13. 44 these 
integrals are split into IR-singular (IR) and IR-finite (fin) parts, 

UM\) = J™ + Jf\ (135) 

The IR-singular parts depend on the scheme adopted to regularize soft and collinear 
singularities. The IR-finite parts are scheme independent and free from soft-collinear 
singularities, but can contain ultraviolet poles. 

Let us start with the two-point functions fll32p . Here only J lo (M|) gives rise to 
IR singularities. This integral is split into 

,2 \ 



IR 



MR 



In 



m 



j 



IR 

lo,DR 



jfxn 



'4vr/i 2 V r(l + e) 



'Att^Y r(l + e) 



M 2 W 



+ 1. 



(136) 



We note that within dimensional regularization the UV and IR singularities can- 
cel each other and the massless two-point function vanishes, J{a DR + Jf° = 0. 



The three-point functions Jgfe(M^) and Jnb{M^) are free from IR singularities and 
the singularities originating from J 8 (M|) and Ji (M|) do not need to be consid- 
ered since the coefficients associated with these scalar integrals are of order M\ 
(see Appendix |E]). The remaining three-point functions in (11331) contain soft and 
collinear singularities. For them we find 



J. 



IR 

7,MR 



IK 



DR 



V s 



En 




(137) 



and 



Ji 



IK 



..MR 



Ji 



IR 

9a, DR 
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J 9a 




(138) 
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where Li 2 (x) = — f$ dtln(l —t)/t. The finite and singular parts for Jn (M|) = 
J$a{Mjd\a-+i are constructed in the same way. 

The singular parts of the subtracted four-point functions (j!34p can be related to 
the ones of the three-point functions, 



tIR 

J 12 — 


tIR 

J 13 




tIR 

- J 7 ' 
U 




tIR 


tIR 
J 14fe 


t*->u 


1 tIR 

t 9<I ' 


(139) 



in both regularization schemes. This implicitly defines the remainders as 
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(140) 



Using the explicit analytic expressions for the infrared singular four-point and three- 
point functions [HI 02] we obtain 
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E Explicit result for the virtual corrections 

In this appendix we present explicit analytic expression for the functions H\(M V ) 
defined in (1741) . These functions describe the contribution of the unrenormalized 
Feynman diagrams of Fig. [2] to the unpolarized cross section. They consist of linear 
combinations of the scalar integrals defined in Appendix [Cj 



H{(M V ) = Y.K){M 2 v )Re[J J {M 2 v ] 

3 

The coefficients of the function H^(M V ) read 

K^(My) - 



4S 2 + 3(t 2 + ^) Jj_ 1 



for I = A, N, X, Y. 
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The only difference between H±(My) and the equally named function in Ref. [9J is 
due to the fact that H±(My) in Ref. [9] includes the contribution of the fermionic 
wave-function renormalization constants, which modify the coefficients Kq and K A a 
(see eqs. (54) and (55) in Ref. [Hj). 

For the coefficients of the function Hf(Ml) we obtain 
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- s(t + u) 
KUM 2 V ) = Kl a {M 2 v ) 



(144) 



For My = Mw the function Hf(My) is identical to the equally named function in 
Ref. 0. 

The only non- vanishing coefficients of the function H^(My) read 
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F Real corrections 

In Tableland Table [2] we list the dipoles that were used to calculate the subtraction 
terms in (1911) for the massive regularization and in f 1 1 U 9 j) for the dimensional reg- 
ularization, respectively. We give references to the explicit formulae for the dipole 
terms and the phase-space mappings in the original paper [T7] and [TBI [19] . 
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Dipole 


Type (emitter, spectator) 


eq. no. 




„sub 
yab,r 


massless IS, massless IS 


(3.22) 


(3.25)-(3.27) 


„sub 
UaW,T 


massless IS, massive FS 


(A.l) 


(3.12) 


„sub 
y\Va,r 


massive FS, massless IS 


(A.l) 


(3.12) 


sub 
j ak,T 


massless IS, massless FS 


(3.9) 


(3.12) 


„sub 
Hka,T 


massless FS, massless IS 


(3.9) 


(3.12) 


„sub 
ykW,r 


massless FS, massive FS 


(4.4) 


(4.5) 


„sub 
y\Vk,T 


massive FS, massless FS 


(4.4) 


(4.5) 



Table 1: Dipole subtraction terms from Ref. [T7] used to calculate M^ b in (J9~I1) for 
the massive regularization (IS = initial-state, FS = final-state). 



Dipole 


Type (emitter, spectator) 


eq. nos. 


$2,nm 


u QED 


massless IS, massless IS 


(5.136), (5.145) 
in Ref. [IE] 


(5.137), (5.139), 
(5.140) in Ref. [18] 


U W, QED 


massless IS, massive FS 


(5.71), (5.81) 
in Ref. [19] 


(5.73), (5.74) 
in Ref. [19] 


U lW, QED 


massive FS, massless IS 


(5.40), (5.50) 
in Ref. [19] 


(5.42), (5.43) 
in Ref. [19] 


u k, QED 


massless IS, massless FS 


(5.61), (5.65) 
in Ref. [18] 


(5.62)-(5.64) 
in Ref. [18] 


T) a 

u ^k, QED 


massless FS, massless IS 


(5.36), (5.39) 
in Ref. [IE] 


(5.37), (5.38) 
in Ref. [IE] 


^-yk,W, QED 


massless FS, massive FS 


(5.2), (5.16) 
in Ref. [IS] 


(5.3), (5.7), (5.9) 
in Ref. [IS] 


^jW,k, QED 


massive FS, massless FS 


(5.2), (5.16) 
in Ref. [19] 


(5.3), (5.7), (5.9) 
in Ref. [19] 



Table 2: Dipole expressions from Refs. [TH1 US] used to calculate M^ b in (11091) for 
the dimensional regularization. 
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